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ABSTRACT 

Using a nonlocal time-dependent theory of convection, we have calculated the 
linear non-adiabatic oscillations of the Horizontal Branch (HB) stars, with both the 
dynamic and thermodynamic coupling between convection and oscillations carefully 
treated. Turbulent pressure and turbulent viscosity have been included consistently 
in our equations of non-adiabatic pulsation. When the coupling between convection 
and oscillations is ignored, for all models with T e < 7350, the fundamental through 
the second overtone are pulsationally unstable; while for T e < 6200 all the models are 
unstable up to (at least) the 9th overtone. When the coupling between convection and 
oscillations is included, the RR Lyrae instability strip is very well predicted. Within 
the strip the most models are pulsationally unstable only for the fundamental and 
the first few overtones. The turbulent viscosity is an important damping mechanism. 
Being exclusively distinct from the luminous red variables (long period variables), the 
HB stars to the right of the RR strip are pulsationally stable for the fundamental and 
low-order overtones, but become unstable for some oen of the high-order overtones. 
This may provide a valuable clue for the short period, low amplitude red variables 
found outside the red edge of the RR strip on the H-R diagram of globular clusters. 
Moreover, we present a new radiation modulated excitation mechanism functioning in 
a zone of radiation flux gradient. The effects of nonlocal convection and the dynamic 
coupling between convection and oscillations are discussed. The spatial oscillations of 
the thermal variables in the pulsational calculations have been effectively suppressed. 

Subject headings: Convection - stars: horizontal branch - stars: Oscillation - stars: 
variables; RR Lyrae 

1. Introduction 

The main goal of the series of our present work is to study the effects of convection on the 
stability of variable stars of different types. Gough (1977) and Xiong(1977) have pointed out 
independently that a nonlocal time-dependent theory of convection is needed in order to account 
for the dynamic coupling between convection and oscillations. Based on the general equations of 
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fluid dynamics and the theory of turbulence, we have developed such a nonlocal time-dependent 
statistical theory of convection (Xiong, Cheng & Deng 1997), which gives an equation for turbulent 
viscosity that is formally similar to the Stokes equation for viscous gas. The theoretical reasonings 
and the correct expression of the gradient approximation for the triple correlations (the 3rd-order 
correlation functions of the turbulent fluctuations of velocity and temperature) are provided by 
the dynamic equations. We have presented results of the linear nonadiabatic oscillations for long 
period variables in the first paper of this series (Xiong, Deng & Cheng, 1997, hereafter Paper I), 
in which a so-called Mira instability region outside the Cepheid instability strip has been defined. 
In this paper, we will apply our theory to the study of the pulsational stability of HB stars. 

The theoretical progress of RR Lyrae star research is due to many authors, among which we 
recall Christy (1966), Iben (1971), Baker and Gough (1979), Deupree (1977a,b,c), Xiong (1980), 
Stellingwerf (1984), Bono and Stellingwerf (1993), Guzik and Cox (1993) and Bono et al (1995, 
1997). Christy(1966) and Iben(1971) completely ignored the coupling between convection and 
oscillations, therefore they gave no definition for the red edge of the RR instability strip. Baker & 
Gough (1979) followed Gough's local time-dependent mixing length theory of convection (Gough 
1977), while Xiong (1980) used a local time-dependent statistical theory of convection (Xiong 
1977). The later two papers gave very close results, and both reached a definite red edge of the 
RR instability strip. The detailed comparison of the two was given by Balmforth & Gough (1987). 
The coupling between convection and oscillations was attempted by Stellingwerf (1984) and Bono 
et al (1993, 1995, 1997) using a rather simplified nonlocal time-dependent theory of convection. In 
this paper, we have put great emphasis on the theoretical baseline of the convection theory, using 
which we have calculated the linear non-adiabatic pulsation for RR Lyrae stars. We are going to 
show the advantages of the present theory over the local time-dependent theory of convection. 

It is claimed that a bunch of short period, small amplitude variables were found outside the 
RR instability strip on the H-R diagram of some globular clusters (Yao 1981, 1986, 1987; Yao et 
al. 1993a, b, 1994). In this paper we try to give a theoretical interpretation of such variables. In 
the following section, the main results of the calculated models for non-adiabatic pulsation are 
presented. We also proposed a new excitation mechanism, i.e. the so called "radiation modulated 
excitation mechanism", which works in the region of radiation flux gradient. We have explained 
in Section || the reason for the spatial oscillations of the thermodynamic quantities that happen in 
the numerical computations of stellar pulsation . A concise discussion and conclusions of our work 
are in Section ^| 

2. The Numerical Results 

We have calculated the fundamental through 11th order overtone modes of linear non-adiabatic 
oscillations for six series of HB star models with a intermediate metal abundance. All the model 
series have the same chemical composition: (X,Z) = (0.76,0.0004). Their mass, luminosity, 
chemical composition and the characteristics of pulsational instablity region are listed in Table 1. 
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The model parameters are: M = (0.65 -0.80)M Q ; L = (1.5 - 2.2) x 10 35 erg/s (M bol = 0.77-0.35); 
and effective temperature in the range of 5150 - 80751T, where the coolest model is very close 
to the Hayashi track; the convection parameters c\ = ci = 0.75 — 1.00, which correspond to the 
same energy transport efficiency as taking a = l/H p ~ 1.5 — 2.0 in the original Vitense theory, 
and the e-folding length of convective overshooting of about (1.05 - lA)H p . The coupling between 
convection and oscillations is carefully treated by using our nonlocal time-dependent theory of 
convection (Xiong, 1989, Xiong, Cheng & Deng 1997). All the mathematical and numerical 
schemes are the same as in Paper I to which we refer for all the details in the equations and 
boundary conditions. 

The lower boundary is set in the convective overshooting zone (T = 8 x 10 6 -fT), and the 
upper one is set optically thin enough in the stellar atmosphere (optical depth r = 0.01). We use 
a simplified MHD equation of state (Hummer and Mihalas, 1988; Mihalas et al. 1988; Dappen 
et al. 1988) in the present work. An analytic approach to the OPAL tabular opacities (Rogers 
and Iglesias 1992) and the low-temperature tabular opacities (Alexander 1975) is used for the 
calculation of opacity. The Henyey method (Henyey et al. 1964) is used for numerical calculations. 
The envelope model of stars is divided into 800 zones with approximately equal intervals in In T. 
However, the mesh points have been increased in the stellar atmosphere, where the temperature 
gradient is very small. The variations of the temperature and pressure within any mesh zone are 
less than 0.03 dex. These mesh points have enough spatial resolution for the highest overtone in 
the present work. The eigen-functions of the 10th overtone for a red HB model are plotted in 
Fig. 

2.1. Radiative Modulation Excitation 

In order to demonstrate the effects of convection on the stability of stellar pulsation, we have 
considered separately two cases, i.e. including and ignoring the coupling between convection and 
oscillations (however, convection has been included in the calculation of equilibrium models). The 
treatment of decoupling can be referred to Paper 1. Simply saying, it is by setting the convective 
variables to zero. Table 2 gives the linear amplitude growth rates for the first 12 eigen-modes for 
the case that the coupling between convection and oscillations is ignored. The first column is the 
serial number of the model, the second is the effective temperature T e , 3rd - 14th columns are 
respectively the amplitude growth rates for the fundamental to the eleventh overtone. As it is 
clear from Table 2, for all models with T e < 7325, the fundamental through the second overtone 
are pulsationally unstable; while all the models with T e < 6200 K are unstable up to (at least) 
the 9th overtone. This makes it clear that the red edge of the RR instability strip cannot be 
determined when the coupling between convection and oscillations is ignored. Fig. ^a-f depicts the 
variation of the integrated work versus the depth (log T) for six HB stellar models with different 
effective temperatures. 

It can be seen from Fig. || that the excitation of pulsationmainly comes mainly from 3 regions, 
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namely, the k mechanism (KM) which works in the vicinity of the helium and hydrogen ionization 
regions, and the modulation excitation mechanism working in the gradient zones of radiation flux 
(RGZ) on the top and at the bottom of convection zone, about which we are going to talk in details 
below. For the HB stars with T e < 5800-fT, the ionization region of helium is almost completely 
convective and the radiation energy transport is nearly negligible (L r /L < 1%). Therefore, the KM 
working near the hydrogen and helium ionization regions is not at all important. The pulsation is 
excited by the KM and the radiation modulated excitation (RME) mechanism above the hydrogen 
ionization region, where dW p /dlogT < (Fig. §a,b). There is a lot of work concentrated on the 
KM in the literature. We will focus on the new mechanism of pulsational excitation, i.e. the 
modulation excitation that works in the RGZ on the top and at the bottom of the convective zone. 
In such RGZ the radiation flux will be modulated when the star pulsates. The function of such 
modulation is to gain energy from the surroundings through radiation while the fluid element is 
hot, and to loss energy while it is cool. Thus on the P-V work plot, a positive Carnot cycle forms, 
that is to say, radiation energy is converted into pulsational kinetic energy by such a process. 
In the following, we are going to explain how such an excitation mechanism works by using the 
integrated work. Following our previous work (Cheng & Xiong 1997), the normalized integrated 
work can be expressed as 

7r r M ° 1 

W P = -— / -=/ m (SPSp*) dM r , (1) 
where, E k is the total kinetic energy of oscillations, 

E k = -oj 2 t / Sr5r*dM r , (2) 

2 JM b 

Starting from the equation of conservation of energy, one can easily have, 

w - T ^t t = - {T3 - 1) ''^ (L ' +L ' +L ' h (3) 

where L r , L c and Lt are respectively the corresponding luminosities of the radiative, convective 
and turbulent kinetic energy fluxes. Inside the equilibrium model of stellar envelope, 

L r + L c + L t = L = const. (4) 

Linearizing Eq. we have 

SP = T x -Sp - (r 3 - 1) -At^- (SL r + 5L C + 5L t ) , (5) 
p ilo aM r 
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Inserting Eq. (||) into Eq. (|l|), we have 
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For not losing any generality and for the sake of convenience in discussing the coupling 
between convection and oscillations, we have kept in Eqs. (||) and the terms related with 
convection, i.e. 5L C and 5Lt- If 5L C and 5Lt are set to zero, this will lead to exclusion of the 
coupling between convection and oscillations. We note 
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where a and 5 are the coefficients of expansion and compression of gas. Putting Eqs. (0) - (11) 
into Eq. ©, Wp can be expressed as a sum of the following five terms 



W P = W K m + Wrme + W C e + W C me + Wi 4 , 
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where Wk m and Wrm e are the contributions of the nomal KM and RME respectively, 
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Wce and Wqme are the contributions of the convective thermal flux and the convective (flux) 
modulation excitation (CME), 
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and Wi,t is the contribution of the turbulent kinetic energy flux, 

In our nonlocal convection theory \L\/L\ <C 1, so Wt, t is negligible in comparision with Wkm-, 
Wrme, W C e and W C me- 

When d{[x T + X P ^2/(T 2 - 1) - 4] \5T/T\} /dr > (< 0), KM plays a destability (stability) 
effect for stellar oscillations. The destability effect is enhanced in the hydrogen and helium 
ionization regions, because T 2 is small due to the ionization of hydrogen and helium (Gamma 
mechanism). When the coupling between convection and oscillations is ignored, Wqe = Wqme = 
and Wp = Wkm + Wrme- This theoretical prediction has been justified by our numerical 
calculations. 

The integrated works Wkm and Wrme are also plotted in Fig. |2[ It can be seen from Fig. [2] 
that in the stellar interior, where the pulsational amplitudes increase outwards, the Wkm increases 
(decreases) with increase (decrease) of kt- 

There is a radiative damping region (dWp/dLogT > 0) in the deep radiative interior, where 
Wkm < 0. For stars with very low effective temperature, such a radiative region goes even deeper 
into the stellar interior where the pulsational amplitude is small, and the radiative damping is small 
as well, as this is seen in Fig. ^a,b. As the temperature increases, the envelope convective zone 
retreats towards the surface. The radiative damping region progresses outwards, and the radiative 
damping increases (Fig. §c-f). There is also a radiative damping region [d?W km I dLogT > 0) in 
the stellar atmosphere, because the pulsational amplitude of 5T/T decreases towards surface. 

The RME exists only in the RGZ. In a stellar envelope RGZ is present only at the top and 
the bottom of the convective zone [Up Radiative Gradient Zone (URGZ) and Bottom Radiative 
Gradient Zone (BRGZ)]. In the BRGZ where dL r /dr < and Xt + X P r 2 / (r 2 - 1) - 4 < 0, 
Wrme is positive. Therefore, the BRGZ is an excitation region of oscillations (Fig. Qb-e). In the 
URGZ the above two inequalities are both reversed. This is another excitation region (Fig. |2|a-d). 
Such an excitation mechanism is related to opacity k, but it is distinct from the normal KM. The 
normal KM can work in purely radiative stellar configuration, for example, such as the excitation 
of (3 Cephei variables, which have almost no convection zone. This new exciting mechanism, 
however, works only in the RGZ. We will call it the radiation modulated excitation (RME), so as 
to distinguish it from the normal KM. 

The location of the RME at the top of a convective zone overlaps with that of the normal KM, 
so it is difficult to distinguish them in the total power diagram (Wp — logT). But, they are clearly 
distinguished with the curves of Wkm and Wrme- The location and the efficiency of the RME 
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at the bottom of the convective zone depend very sensitively on the effective temperature of the 
stellar model under consideration. For extremely low temperature models, the envelope convective 
zone is very extended, and this leads to a very deep modulated excitation region at the bottom of 
the convective zone. However, the pulsational amplitude in this region is very small, so a very 
weak excitation is expected. The main excitation comes from the KM and RME mechanisms in 
the URGZ (Fig. §a,b) for the cool HB stars. In Fig. ||b, the bottom RME is only barely visible. As 
the effective temperature increases, the bottom RME raises up together with the shrinking of the 
convective zone (Fig. ||a-d). The effect of bottom RME is clearly illustrated in Fig |2|b and Fig ^c. 
The bottom RME starts to work at logT ~ 5.7 and 5.4 respectively, corresponding to the location 
of the bottom of the convective zone and being much lower than the second ionization region of 
helium where the normal KM works. For a HB star with T e = 5900 the bottom RME region is 
connected or overlaps with the second ionization region of helium (Fig. ^c). Numerically, RME is 
stronger than the KM in the second ionization region of helium. For models with T e = 6200JT the 
bottom RME region is near the second ionization region of helium (Fig. ^d). From the total power 
diagram (W P - logT) it is difficult to distinguish the RME from the KM. The effects of RME still 
can be seen from the large increase of integrated work Wrme near logT = 4.8 corresponding to 
the location of the bottom of the convective zone. 

For the blue HB stars , the convective energy transport is insignificant since the convective 
envelope is very shallow, and the KM functioning in the second ionization region of helium 
becomes dominant, while RME is almost not at all working. We can also see from Fig. ^f that the 
excitation due to the KM cannot balance the radiative damping in deeper layers for high enough 
effective temperature models and the star becomes pulsationally stable. This is the reason of the 
existence of the blue edge of the RR instability strip. 

One can see from Fig. ^ that there is a bump near logT ~ 5.3. It results from the bump of 
the new OPAL opacity. 

2.2. The RR Instability Strip 

As we have repeatedly stressed, when the coupling between convection and oscillations is 
ignored, all the HB star models with T e < 7550-fC are pulsationally unstable. Therefore, no 
definition for the red edge of RR strip can be given. However, when the coupling is taken 
into consideration, the situation becomes very promising. We have listed in Table 3 the linear 
amplitude growth rates of pulsation (per period) for the first twelve modes. From this table we 
find that when the coupling between convection and oscillations is considered, the fundamental 
and low-order overtones (n < 3) become pulsationally stable for all the HB stellar models with 
T e < 5940. Namely, there is an apparent red edge of the RR instability strip. The amplitude 
growth rates of the fundamental and the first overtone versus the effective temperatures for the 
model series 2-4 are plotted in Fig. ||. In Fig. ||| the integrated works versus depth are given 
for the fundamental modes of 6 HB stellar models, where Wp, Wpt and W v i s are, respectively, 
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the contributions due to the gas (including radiation), turbulent pressure and turbulent viscosity, 
while W^i = Wp + Wpt + W v i s is the total integrated work, the numerical value of which at the 
surface is approximately the same as the linear pulsational amplitude growth rate r] = —2TTL0i/L0 r 
computed with the non-adiabatic scheme. Here uji and uj t are the imaginary and real part of 
the complex pulsational angular frequency (a; = iooi + u r ). For the real definitions of Wp, Wpt 
and Wvis we refer to our previous work (Cheng & Xiong 1997). Here we only recall that of W v i s 
to make our discussion clearer. — W v { s is the random turbulent kinetic energy converted from 
the pulsational energy due to turbulent viscosity, which is finally turned into thermal energy 
through molecular viscous processes. W v i s is always negative. Therefore, This implies a damping 
mechanism. 

We can see from Fig. || that, for HB models with T e > 6000-RT, the radiation energy 
transport is becoming dominant gradually in their second ionization region of helium. For 
these stars, the main source of excitation is the KM in the second ionization region. When the 
effective temperature drops below 5900K for the HB models, the hydrogen and helium ionization 
regions are fully convective, and the radiation contributes very little to the total energy flux 
{L r /L < 1%). Therefore, the KM is almost inhibited. The convection becomes the dominant 
factor for determination of the pulsational stability of stars. Thus, the low temperature HB stars 
with T e < 5940K become pulsationally stable. The coupling between convection and oscillations is 
the cause why there exists a red edge for the RR Lyrae instability strip. Fig |6| plots the variations 
of temperature for the lower boundaries of the convective zone and the ionization regions of 
hydrogen and helium versus the effective temperature of the HB models under consideration. As 
demonstrated in Fig. ||, the red edge of the RR Lyrae instability strip is defined when the bottom 
of the convective zone is just extended slightly below the second ionization region of helium. 



2.3. The Short Period Red Variables Beyond the RR Lyrae Strip 

Yao and his collaborators found some short period red variables outside the RR Lyrae strip 
in the H-R diagrams of a few globular clusters (Yao 1981, 1986, 1987; Yao et al 1993a,b, 1994). 
The light amplitude of such red variables is as small as a few percent of a magnitude. These stars 
are usually multi-periodic or irregular; their typical period or the time-scale of variability is a 
few tenths of a day. Taking the variation of light of such red variables as a result of intrinsical 
pulsation, the natural conclusion is that the stars all pulsate at high-order overtones. In general, 
the dwarfs tend to be non-radially excited while the giants are radially excited. Therefore, we think 
that the short period red variables are pulsating in high radial overtones. Of course, we cannot 
exclude the possibility that they pulsate at low degree non-radial modes. It is very interesting to 
see from Table 3 that, when the coupling between convection and oscillations is considered, all 
the fundamental and low-order overtones of the red HB models for T e < 5940-fT are pulsationally 
stable; while some high-order overtones (n > 4)are pulsationally unstable. We also show that 
these low temperature red variables possess nearly the same period as that of the variables found 
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to right side of the RR Lyrae strip by Yao and his collaborators. A general regularity of variable 
stars is that when a star pulsates at high-order overtones, its pulsation amplitude is usually 
very small and tends to be multi-periodic. The results of our theoretical modeling are likely to 
be a reasonable explanation for the red variables found outside the RR Lyrae strip in globular 
clusters and for their qualitative properties, though there still exist some questions to be studied. 
According to our theory, it seems that all the low temperature HB stars ought to be pulsationally 
unstable at high-order overtones. This raises a naive question: why these short period oscillations 
were observed only for some, not for all the low temperature HB stars? A straightforward answer 
to such a question is that they all might be pulsationally unstable but are not detectable with 
current instrumental precision. This seems to be a reasonable answer. However, there is still a 
question to be asked: what is the drive to make some of the stars pulsate with higher amplitudes? 
We are not ready to give satisfactory explanations to such questions. 



2.4. Comparison with Other Theoretical Results and Observations 

The RR Lyrae instability strip yielded by our numerical computations of the linear non- 
adiabatic oscillations is 3.774 < logT e < 3.884 for the second model series {c\ = C2 = 0.85). The 
location of instability strip of the fundamental mode agrees well with Bono's nonlinear theoretical 
results (Bono, et al. 1997). The temperature of the red edge for the fundamental mode is slightly 
higher, while the temperature of the blue edge is slightly lower than theirs. Our blue edge of 
the first overtone is much hotter than that given by Bone's nonliner results. Their blue edge, we 
think, seems to be too red . Our theoretical red edge of the first overtone is slightly hotter than 
the fundamental mode red edge. Obviously, it is too red in comparison with observations. It is the 
intrinsic defect of the linear theory. This is a problem of mode choice. It is impossible to decide 
which mode can finally become an unstable one with limited amplitude for the linear theory when 
several modes are linearly pulsationally unstable. 

For the globular cluster M3 whose metallicity is close to what we adopted in this work. 
Sandage (1990) determined observationally a RR Lyrae instability strip 3.791 < logT e < 3.880. 
This infers that the theoretical result matches very well the observed one for the blue edge of 
the instability strip, while the red edge is 240 K lower than the observed one. For the same 
mass, luminosity, effective temeperatute and chemical composition, the convection zone deepens 
as ci increases. Therefore, the location of the red edge of the instability strip can be adjusted by 
modifying the convective parameter c\ . We calculated the linear non-adiabatic oscillations of three 
series of HB star models with different c\. The effective temperatures of the blue and red edges 
for the fundamental and the first overtone are listed in Table 1. The red edge is still 170 K lower 
than observation even when c\ is increased to 1.0. The energy transport efficiency in this case 
corresponds to that of taking a = l/H p ~ 2.0 in the original Vitense theory. It seems impossible 
that ci can be much larger than 1. We do not know whether this is due to the imperfection of our 
theory or because the temperature of the red edge determined by Sandage is too high. 
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2.5. Luminous Red Variables (LRV) and Red HB Stars (RHBS) 

Although LRV and RHBS are in the different evolutional statuses, they all possess very 
extended convective envelopes. Convection is the most dominant factor controlling their pulsational 
stability However, their pulsational properties are entirely different: LRV are pulsating at the 
fundamental and the first few overtones, while all higher (n > 4) overtones are stable. On the 
contrary, for RHBS the fundamental and lower order overtones (n ~ 4) are pulsationally stable, 
while some higher order overtones are unstable. It is hard to understand this acute difference at 
the first glance. In fact, this is because these two types of stars have quite different mass-luminosity 
ratios: The masses of LRV do not differ much from those of RHBS, while the luminosities of LRV 
are higher than those of RHBS up to 2 orders of magnitude. Thus their mass-luminosity ratios 
differ about in 2 orders of magnitude. This implies that the envelope structure of LRV is very 
different from that of RHBS. For LRV, the central condensation is very high. This will lead to 
very different pulsational properties for these two types: 

1. An important factor is that the thermal time scales of pulsational amplitude growth for these 
two types are very different, and this is due to the tremendous difference of mass-luminosity 
ratios. According to Eq. (Q) in Paper I, this thermal time scale can be expressed as 
Airr^pCpT/L. Since the luminosity of LRV is high and the gaseous density in their envelope 
is very low in comparison with RHBS, the thermal time scale of LRV is very short (compared 
with their pulsation period). It has been shown by actual calculations that the amplitude 
growth rate (the amplitude increment in a pulsational period) of LRV is higher than that 
of RHBS by about 2-3 orders of magnitude. This means that, for LRV, the excitation 
(damping) by radiative (and convective) energy transport is very strong. Therefore, for 
fundamental and low-order overtones, the effect of turbulent viscosity will be relatively 
weak (please refer to Fig. |3] in Paper I). Besides, LRV will tend to pulsate at fundamental 
and low-order overtones because lut c <C 1 in most regions of LRV interiors. Here r c is the 
convective inertial time scale. As the frequency increases, the effects of turbulent viscosity 
increase rapidly [the turbulent viscosity is proportional to u>r c /(l + lo 2 t 2 )]. Thus their high 
overtones tend to be pulsationally stable (prease refer to Fig. ||| in Paper I). For RHBS, since 
the amplitude growth rates are about 2-3 orders of magnitude less than those of LRV, the 
effects of turbulent viscosity become relatively important (see Fig. |J), and will damp the 
pulsation. The tremendous difference of mass-luminosity ratios is another reason for the 
distinct properties of the two types of low temperature variables. 

2. The distributions of lut c are totally different in the interiors of these two types. In fact, the 
interaction between convection and oscillations is very complicated. The effects of convection 
on stellar pulsation depend not only on the variation amplitudes of convective variables, but 
also on their relative phases. Both the amplitude and phase depend critically on the ratio of 
convective inertial time scale (r c ) and that of pulsation, i.e. u>t c . Comparing the integrated 
work plot of the present work (Fig. [|) with that of long period variables in our previous paper 
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(Paper I), we can see that, the quantity ut c is very different in the interiors of these two 
types of stars. For the LRV, lot c < 1 in most part of the stellar interior except the very outer 
part near stellar surface; while for the RHBS ujt c 3> 1 in almost the whole stellar interior. 
Referring to Fig 3 of our previous work (Paper I), it is clear that the hydrogen and helium 
first ionization regions are the main excitation region of stellar pulsation for the LRV, where 
convection tends to destabilize pulsation, while the turbulent viscosity is the main damping 
mechanism for high-order overtones, due to turbulent viscous dissipation (oc wr c /(l + uj 2 t 2 )). 
When ujT c <C 1, the turbulent viscosity is proportional to ujt c , which increases very quickly 
with increasing pulsational frequency. It becomes opposite for the RHBS, for which ujt c 3> 1 
holds in most part of the stellar interior. The excitation due to the first ionization region 
of helium and a part of ionization region of hydrogen is restrained, and the damping effect 
of turbulent viscosity on high-order overtones also diminishes. This tends to stabilize the 
fundamental and the low-order overtones of the RHBS while the high-order overtones become 
pulsationally unstable. Following such discussion we conclude that the difference in lot c is a 
reason for the distinct properties of the two types of low temperature variables. As for the 
question of why cut c 's are so different in the two type of stars, we understand that it is due 
to the convective inertial time scale r c oc Hp/x. Here Hp = r 2 P/GM r p cx r 2 P/M r p is the 
local pressure scale height and x is the r.m.s. component of the turbulent velocity. Relative 
to the RHBS, the LRV have higher central condensation. Towards stellar center, r goes 
much faster with increasing temperature. Fig. [6| plots r, Hp and wr c of fundamental mode 
versus temperature for a RHBS model and a LRV model respectively. It is clearly shown 
in Fig. fj] that there exists very large difference between the two types of low temperature 
variables of different luminosities. Following the above discussion, we can easily understand 
the reason of the distinct properties of the two types of low temperature variables although 
their variabilities are both affected critically by convection. Fig. |5| plots the variations of the 
integrated work of the first, second, and fourth overtones versus depth for a T e = 5150i"C 
RHBS model. From this we see that, unlike luminous red stars, the turbulent viscosity 
component W v i s does not increase quickly with uj r . 

3. The Spatial Oscillations of the Convective Variables 

In the theoretical calculation of the steller oscillations, if the local theory is used to treat 
the coupling between the convection and pulsation, the troublesome spatial oscillations of 
thermodynamic variables will occur (Keeeley, 1977; Baker and Gough, 1979; Gonczi and Osaki, 
1980). These spatial oscillations are due to the local treatment of the convection, so we could 
expect that it will automaticly disappear if the non-local treatment of convection is adopted. 
Unfortunately, there still exist intense spatial oscillations of convective and thermodynamic 
variables in our previous work on the solar radial p-mode oscillations, although a non-local 
convection theory was used (Cheng and Xiong, 1997). Now, we use an asymptotic analysis to 
inspect the cause of spatial oscillations. It is necessary to review briefly the main points of our 
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non-local convection theory first. Our convection theory is a statistical theory of correlation 
functions. We derived a set of dynamic equations of the auto- and cross-correlations of the 
turbulent velocity and temperature fluctuations. For example, the dynamic equation for the 
auto-correlations of the turbulent velocity is as follows (Xiong, Deng and Cheng, 1997), 



__ w H w rj _|_ uu' l w' k X7kui + w'^w^VkW 1 



. ik w'3T' ik w n T\,„ = D« fc . 



pV k (pu' k w H w'i) = -AVZrieXw'iw'i /l e , (18) 



where x = yw' i w h /3 and l e is the average size of the energy-containing eddies. We assume that 
l e ~ c\H p {H p is the local pressure scale height). The right-hand side term of Eq. (|l8|) is the 
turbulent dissipation, which determines the efficiency of the convective energy transport. Thus, c\, 
which is a convective parameter relevant to the turbulent dissipation, functions like the parameter 
a = l/Hp in the mixing-length theory (Bohm-Vitense, Z. 1958). In our theory, the efficiency of the 
convective energy transport is approximately the same as letting a = l/H p = 2c\ in the oringinal 
mixing-length theory. 



pw h w' J is the turbulent Reynold stress. We can decompose it into an isotropic component 
(turbulent pressure) and a nonisotropic one (turbulent viscosity). An expression for the turbulent 
viscosity very similar in form to the formula of Stokes viscosity has been derived (Xiong, Deng 
and Cheng, 1997). To reveal the nature of the non-local convection, we consider only the isotropic 
component of Eq. (|l8|), i.e. 



3Dx 2 



-x 2 V k u k + aV k ( BUk 



2 Bt v Dt 

+ pVkipu'kw'tw'*) - 2V3r] e x 3 /l 



(19) 



where V k = w' k T' jbT embodies the correlation of the turbulent velocity and temperature 
fluctuations. Eq. ( |l9| ) is the dynamic equation of turbulent pressure. But it can also be considered 
as the equation of conservation of turbulunt kinetic energy. The left-hand side of Eq. ( |l~9| ) is 
the rate of turbulent kinetic energy, which equals the sum of the terms on the right-hand side. 
The first one is the rate of regular kinetic energy of averagy motion converted to turbulent 
kinetic energy. The second term is the rate of thermal energy converted to turbulent kinetic 
energy through the work done by buoyancy. The last term is the rate of turbulent kinetic energy 
converted to thermal energy due to the molecular viscosity. The triple correlations in the brackets 
of the third term are the turbulent kinetic flux, which represents the non-local turbulent kinetic 
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energy transport by convection. These triple correlations represent all the effects of non-local 
convection. They reveal the essential difference between local and non-local theory: If one ignores 
all the triple correlations, the theory will return to the local convection theory (Xiong, 1977). In 
our earliest works, we used a simplified gradient-type diffusion approximation (Xiong, 1979), i.e. 



vfcu/su)* = -IxVkw'jW 11 = -3/xV fc x 2 , (20) 

where I is the diffusion length. We assume I = c 2 V3H p /4. C2 is another convective parameter 
relevant to the non-local turbulent diffusion. The e-folding length of convective overshooting is 
about lAy/cxc%H p . Let C2 = 0, then we get back to local convection theory. 



Eq. (20) is introduced only as a physically reasonable assumption. Recently, we have 
derived a set of dynamic equations for triple correlations (Xiong, Deng and Cheng, 1997). These 
equations are very complex: The complete set of the radiative fluid dynamic equations are up 
to the 14th order. In the dynamic equations for triple correlations, v! k w' i w h contains not only 
the second correlations of turbulent velocity, but also the cross- correlations of turbulent velocity 
and temperature. For simplicity, omitting the non-corresponding second correlations, we get a 
simplified time-dependent equation for triple correlations: 



ll 



^-i — + u'w'w H = -xlV k w'w H = -3xlV k x 2 . (21) 
Bt 

In our calculation of the solar radial p-modes oscillations (Cheng and Xiong, 1997), a set 



of formulae similar to Eq. (21) had been used to treat all triple correlations (u^T' 2 /T 2 and 
vLw'^T' /T). But the results are just contrary to our expectation: The strong spatial oscillations of 
the convective and related thermal variables existed persistently. 



For radial pulsation, Eqs. (|19|) and (|2l|) can be rewritten as follows 



3Dx 2 



2 Bt 



+ 3 



p Bt 



dM r 



c\r 2 P 



DLi _ Dln(pr 2 ) 
Li 



Bt 



Bt 



t c d ( u r \ Bu r 
Idhxr 7 + ~Bt 



0, 



4v / 37r 2 C2P-Pr 6 x 3 91nx 



GM r 



dM r 



(22) 
(23) 



where 3Li = 2irr 2 pu' 1 w' i w n is the turbulent kinetic energy flux. In the region where ui r r c ^> 1, the 
convective variables x and V change much faster than r, p, T and P of the average fluid fields. 
Therefore, under the first approximation for the purpose of studying the asymptotic properties of 
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the convective variables, we can ignore the slow variations of r, p, P and T of the fluid fields, and 



simply take them as constants. Bearing this in mind and linearizing Eqs. (22) and (|23|), we have 
the following approximated and linearized equations of pulsation, 



\/3c 2 P3: 3 dlnx d f 8L X 
AGM r p dlnrdlnr V Li 



GM r 



V3 Ve (2 + 



+ 



1 6Li 



\^3c 2 r 3 Px 3 dlnx 



lUTc) 



pr 3 L\ dlnr 

px 3 Sx 



C\P x 




(24) 



^/3irc 2 r 3 Px 3 ( d 



GM r 



.din 



r \ x 



Sx\ 
— + 



dlnx 
dlnr 



Jx SL X 

3 (1 + wt c ) —— 

x L\ 



(25) 



Besides the convective variables Sx/x and 5L\/Li, Eq. (24) still contains another convective 
variable 5V/V . Its linearized equation will have other variables, namely Sx/x and 5Z/Z, where 
Z = T' 2 /T 2 represents the auto-correlation of the turbulent temperature fluctuation. In this case, 
the three auto- and self-correlations of the fluctuations of turbulent velocity and temperature are 
closely coupled. We have to solve the 6 linear differential equations for the asymptotic solutions, 
and this is an extremely difficult procedure. Without losing any generality, we will omit the 
terms containing the non-diagonal element 5V/V in Eq. ( |2~4| ) when we discuss Eqs. ( |25| ) and (|2~i|). 
Replacing the variable with 7] = In P, the two equations can be simplified as 



drj 



;i + icoT C ) y 2 ] = o, 



(26) 



dY 2 ( dlnp dlnr\ 4rj e . 

— h 3a H h 2— — Y 2 (2 + iujt c ) Y\ = 0, 

aij V dry arj J c\c 2 ot 



(27) 



where 



dlnx T ^ Sx 5L\ 

a = -T] — pi Yi = —, Y 2 = — 
dlnr x L\ 



(28) 



Neglecting the variation of pr in space and making the variable conversion 



4^2 + ^t c 1 x = ^ 
y c\c 2 1 + iujT c a a 



(29) 
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Eqs. (|2q) and (ETn can be written as 



+ 3ali - (1 + iuJT c ) (2 + iooT c )X = 0, (30) 
dr/ v 

dX 



— + 3aX - /3J(1 + iuoT c ) (2 + wt c )Y x = 0. (31) 
ar\ v 



The characteristic values for Eqs. fl27j) and ( |28| ) are 

A ± = -3a ± pyj (1 + iwr c ) (2 + iujt c ) « -3a ± Q + iwr c ^ /3, (32) 

and their general solutions are 

Yi = aiP" 3o+ l /3 cos (/?cjr c 7? + ^) + a 2 P~ 3a ~% P cos (/?cjr c 7? + 2 ) , (33) 

Y 2 = - \a 1 p- 3a+ ^ cos (^T c r/ + ^) - aaP- 3 ""^ cos (^r c r/ + 2 )1 , (34) 
a L J 

where ai, a 2 , 4>i and 02 are constants. Apparently, the solutions Eqs. ( p3| ) and ( |3"4| ) are of the 
oscillating type, and the wavelength of the spatial oscillation is 

27T -H P « 4.68^. (35) 



(3UJT C UJT C 

Within the convective overshooting zone, it follows from the asymptotic properties of convective 
overshooting (Xiong, 1989) that 



^ / l + 2^7/ e , 0-70175 

a f« ±W = = ± — (36) 

V 3\/3ciC2 

In the upper overshooting zone of a convective region, a is the positive solution. To avoid a too 
fast outwards growing amplitude, we choose A = A+, i.e. we take a 2 = in Eqs. ( |33| ) and fl34|). 
On the other hand, in the lower overshooting zone of a convective region, where a < 0, we take 
A = A_, i.e. ai = 0. 



According to Eqs. (p^), (p3|), (|34"D and (36), the convective variables tend to have an oscillating 



solution with a growing amplitude towards convective overshooting zone, the e-folding length of 
which is 



-3<J ± |/? 



°-^H P . (37) 
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There exist the similar spatial oscillating behavior as above for the auto-correlation of 
temperature fluctuation Z and the cross-correlation V of velocity and temperature fluctuations. 
The procedure to follow is also the same as described above, so we will not go into theirs details. 

Now, it is not difficult to understand why there exist intense spatial oscillations in the region 



where oj t t c 1. Let's look back at Eq. (pip. The triple correlations u^w^w' 1 must approach 
to zero when u) r r c S> 1. That is to say, the nonlocal convection theory approaches to the local 
one. Thus the spatial oscillations occurring in the calculation of local convection theory will 
certainly come out. To avoid spatial oscillations, we go back to the simple gradient type diffusion 
approximation respectively by Eq. (|20|). The corresponding linear pulsational equation becomes 

+ 3ali - aY 2 = 0. (38) 

drj 

Making a new variable conversion 



Y 2 = J^(2 + iuT c )-X=^X, (39) 
V c\C2 a a 



Eqs. (|38j) and ([271) can be written as 

+ 3aFi - f3W2 + iuJT c X = 0, (40) 

drj 
dX 



— + 3aX - (3W2 + iwT c Yi = 0. (41) 
drj 

Their characteristic values are 

A± = 3a ± V2 + iuJT c (5i = 3a ± \Z2f3 x (cos (j) + i sin , (42) 

where 2<p = arctan^^. It is clear that if ujt c — ► oo, we have 4> — > ir/4. Therefore, in the 
upper overshooting zone where a > 0, we chose the solution with A = A^; while in the 
bottom overshooting zone where a < we chose the solution with A = X'_. Then the spatial 
oscllations become solutions attenuating rapidly towards the overshooting zone. By doing so, we 
restrained effectively the spatial oscillations in our calculations. The trace of spatial oscillations 
of eigen-functions can still be found in the outmost zones where u>t c 3> 1 in Fig. [I], and especially 
for convection variables (2/5-2/10)- The eigen- function of turbulent velocity in the outmost stellar 
zones is plotted on Fig. [l]f, where the spatial oscillations are very clear. 

These asymptotic solutions were used as the convective boundary conditions at the top and 
bottom in the numerical calculations of non-adiabatic oscillations. 
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4. Discussion and Conclusions 

We have carried out theoretical calculations of linear non-adiabatic oscillations for six series 
of HB stellar models. The main results can be summarized in the following, 

1. When the coupling between convection and oscillations is ignored, it is not possible to 
interpret the red edge of the RR Lyrae instability strip. All the HB models with T e < 6200i^ 
are pulsationally unstable from the fundamental up to (at least) the 9th overtone; 

2. The intrinsic reason for the existence of the red edge of the RR instability strip is the 
coupling between convection and oscillations. Our numerical results show that, when taking 
the coupling into consideration, all the low temperature HB models with T e < 5940-fT have 
pulsationally stable fundamental and low-order overtones (n < 3). The turbulent viscosity is 
an important damping mechanism of oscillations. 

3. For the low temperature red HB star models some of the high-order overtones (n > 4) 
become pulsationally unstable. The pulsational periods of these high-order overtones are of 
the order of a few tenths of a day. This might be the clue for the red variables outside the 
RR Lyrae instability strip in globular clusters found by Yao and his collabolators(Yao, 1981, 
1986, 1987; Yao et al. 1993a,b, 1994). 

4. The pulsational properties of RHBS are very different from those of LRV. In the current 
works, all low temperature red models with T e < 5940-K" have pulsationally unstable 
fundamental and the low-order overtones (n < 3), while some of their high-order overtones 
(n > 4) are pulsationally unstable. Opposite to this, the LRV outside the Cepheid instability 
strip (i.e. the long period variables) are pulsationally stable for all the high-order overtones 
(n > 4), while the fundamental and the first overtone are possibly unstable. The tremendous 
difference of mass-luminosity ratios is the reason for the distinct properties of the two types 
of low temperature variables. 

5. The spatial oscillations of the convective and the thermal variables in the local time- 
dependent convection theory are effectively restrained in the present paper, where a nonlocal 
time-dependent convection theory is used to treat the coupling between convection and 
oscillations. 

Comparing Tables 2 and 3 we have found that, the effects of convection on the instability 
of such stars still cannot be ignored for the blue HB stars (T e > 7000-fQ, although the 
convective envelope is already very shallow and the convective energy transport becomes negligible 
(L c /L <C 1). These effects are mainly due to the dynamic coupling between convection and 
oscillations, other than the thermodynamic coupling, because in the outermost layer of the blue 
HB stars, a certain amount of turbulent pressure remains, even if the convective energy transport 
is negligible there. We made a numerical test, in which the dynamic coupling was ignored. The 
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results show that all the low-order overtones become pulsationally stable. We are reluctant to 
draw any further conclusion on the pulsational properties of the blue HB stars because there are 
some undetermined factors, i.e. 

1. There are spatial oscillations in the stellar atmosphere (Fig. [l|f). We cannot estimate 
quantitatively their influence on the results of numerical calculations. Their effect, we think, 
may be unimportant to the cool RHBS, because they have extended convective envelopes. 
There are no any spatial oscillations in the main part of the convective zone except the 
atmosphere (Fig. ||). However, convection appears only in a limited surface region and 
the pulsational stability coefficient is small for the blue HB stars. The spatial oscillations 
become, maybe, more important in comparison with the RHBS. 

2. The undetermined structure of the turbulent atmosphere of stars. The atmosphere has an 
unnegligible influence on the pulsational stability of stars, especially for the blue HB stars, 
which have small pulsational stability. It is very difficult to construct an exact enough model 
of turbulent atmosphere. 

This problem deserves future investigation. The high-precision photometry of variables is 
available for us to determine whether there are low-order unstable modes in the blue HB stars. 

The present work is supported in part by the contract 19573018 of the National Natural 
Science Foundation of China(NSFC). 
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Figure Captions 

Fig. 1. — The normalized eigen- functions of the 10th overtone for the coldest HB star model in 
series 2. 

Fig. 2. — The integrated work Wp versus depth for six different HB star models. The coupling 
between convection and oscillations is ignoreed. Wkm and Wrme are respectively the contribution 
components of the k mechanism (KM) and radiative modulation excitation (RME).The radiative 
energy flux L r /L, the adiabatic index T2, and the partial differential of opacity with respect to 
temperature, \ T are also plotted. Three horizontal lines indicate the locations of the ionization 
regions of H and He. 

Fig. 3. — The variation of the growth rates of the fundamental and the first overtone versus 
effective temperature for three model series of HB stars. The solid and long-dashed lines show 
the cases with consideration of the coupling between convection and oscillations (WCC), while the 
dotted and the dot-dashed lines shows the cases ignoring the convection coupling (W - OCC). 

Fig. 4. — The same as Fig. |2[ but the coupling between convection and oscillations is considered. 
Wp, Wpt and W v i s are, respectively, the contributions due to gas (and radiation) pressure, turbulent 
pressure and turbulent viscosity, while W a u is the total integrated work, the other symbols are the 
same as in Fig. ^. 

Fig. 5. — The integrated work virsus depth for the first, second, and fourth overtones of a 
T e = 5150-ftT HB star model, the coupling between convection and oscillations is considered. Refer 
to Fig. H for details. 

Fig. 6. — The bottom temperatures of the hydrogen, the first and the second helium ionization 
regions and convective region versus effective temperature for the model series 2 of HB stars. The 
horizontal lines are the location of the RR Lyrae instability strip. The horizontal line indicates the 
locations of the RR Lyrae instability sttip. 

Fig. 7. — The variation of radius r, mass M r and the local pressure scale height Hp versus depth 
(logT). Panel a.: the envelope model of a HB star with T e = 5300-fC; Panel b.: a long period 
variable envelope model with M = 1M , L = 5000L o , T e = 2700if, x = 0.70 and z = 0.020. 
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